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Non scaling Fixed-Field Alternating Gradient (FFAG) accelerators have an unprecedented po- 
tential for muon acceleration, as well as for medical purposes based on carbon and proton hadron 
therapy. They also represent a possible active element for an Accelerator Driven Subcritical Reactor 
(ADSR). Starting from first principle the Hamiltonian formalism for the description of the dynamics 
of particles in non scaling FFAG machines has been developed. The stationary reference (closed) 
orbit has been found within the Hamiltonian framework. The dependence of the path length on the 
energy deviation has been described in terms of higher order dispersion functions. The latter have 
been used subsequently to specify the longitudinal part of the Hamiltonian. It has been shown that 
higher order phase slip coefficients should be taken into account to adequately describe the acceler- 
ation in non scaling FFAG accelerators. A complete theory of the fast (serpentine) acceleration in 
non scaling FFAGs has been developed. An example of the theory is presented for the parameters 
of the Electron Machine with Many Applications (EMMA), a prototype electron non scaling FFAG 
to be hosted at Daresbury Laboratory. 

PACS numbers; 29.20.-c, 29.20.D-, 41.85.-p 
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I. INTRODUCTION 

Fixed-Field Alternating Gradient (FFAG) accelerators 
were proposed half century ago [U, H, 0] , when accelera- 
tion of electrons was first demonstrated. These machines, 
which were intensively studied in the 1950s and 1960s but 
never progressed beyond the model stage, have in recent 
years become the focus of renewed attention. Accelera- 
tion of protons has been recently achieved at the KEK 
Proof-of-Principlc (PoP) proton FFAG [|. 

To avoid the slow crossing of betatron resonances as- 
sociated with a typical low energy-gain per turn, the first 
FFAGs designed and constructed so far have been based 
on the " scaling" principle. The latter implies that the or- 
bit shape and betatron tunes must be kept fixed during 
the acceleration process. Thus, magnets must be built 
with constant field index, while in the case of spiral-sector 
designs the spiral angle must be constant as well. Ma- 
chines of this type use conventional magnets with the 
bending and focusing field being kept constant during 
acceleration. The latter alternate in sign, providing a 
more compact radial extension and consequently smaller 
aperture as compared to the AVF cyclotrons. The ring 
essentially consists of a sequence of short cells with very 
large periodicity. 

Non scaling FFAG machines have until recently been 
considered as an alternative. The bending and the focus- 
ing is provided simultaneously by focusing and defocus- 
ing quadrupole magnets repeating in an alternating se- 
quence. There is a number of advantages of the non scal- 
ing FFAG lattice as compared to the scaling one, among 
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which are the relatively small transverse magnet aperture 
(tending to be much smaller than the one for scaling ma- 
chines) and the lower field strength. Unfortunately this 
lattice leads to a large betatron tune variation across the 
required energy range for acceleration as opposed to the 
scaling lattice. As a consequence several resonances are 
crossed during the acceleration cycle, some of them non- 
linear created by the magnetic field imperfections, as well 
as half-integer and integer ones. A possible bypass to this 
problem is the rapid acceleration (of utmost importance 
for muons), which allows betatron resonances no time to 
essentially damage beam quality. 

Because non scaling FFAG accelerators have other- 
wise very desirable features, it is important to investi- 
gate analytically and numerically some of the peculiari- 
ties of the beam dynamics, the new type of fast accelera- 
tion regime (so-called serpentine acceleration) and the ef- 
fects of crossing of linear as well as nonlinear resonances. 
Moreover, it is important to examine the most favorable 
phase at which the cavities need to be set for the optimal 
acceleration. Some of these problems will be discussed in 
the present paper. 

An example of the theory developed here is presented 
for the parameters of the Electron Machine with Many 
Applications (EMMA) 0] , a prototype electron non scal- 
ing FFAG to be hosted at Daresbury Laboratory. The 
Accelerators and Lasers In Combined Experiments (AL- 
ICE) accelerator 0] is used as an injector to the EMMA 
ring. The energy delivered by this injector can vary from 
a 10 to 20 MeV single bunch train with a bunch charge of 
16 to 32 pC at a rate of 1 to 20 Hz. ALICE is presently 
designed to deliver bunches which are around 4 ps and 
8.35 MeV from the exit of the booster of its injector line. 
These are then accelerated to 10 or 20 MeV in the main 
ALICE linac after which they are sent to the EMMA 
injection line. The EMMA injection line ends with a 
septum for injection into the EMMA ring itself followed 
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by two kickers so as to direct the beam onto the correct, 
energy dependent, trajectory. After circulation in the 
EMMA ring, the electron bunches are extracted using 
what is almost a mirror image of the injection setup with 
two kickers followed by an extraction septum. The beam 
is then transported to a diagnostic line whose purpose it 
is to analyze in as much detail as possible the effect the 
non scaling FFAG has had on the bunch. 

The paper is organized as follows. Firstly, we review 
some generalities and first principles of the Hamiltonian 
formalism [S] suitably modified to cover the case of a non 
scaling FFAG lattice. Subsequently the synchrobetatron 
framework is applied to determine the energy dependent 
reference orbit. Stability of motion about the stationary 
reference orbit is described in terms of betatron oscilla- 
tions with energy dependent Twiss parameters and be- 
tatron tunes. Dispersion, measuring the effect of energy 
variation on the path length along the reference orbit is 
an essential feature of non scaling FFAGs. Within the 



developed synchrobetatron formalism higher order dis- 
persion functions have been introduced and their con- 
tribution to the longitudinal dynamics has been further 
analyzed. Finally, a complete description of the so-called 
serpentine acceleration in non scaling lepton FFAGs is 
given together with conclusions. The calculations of the 
reference orbit and phase stability are detailed in the ap- 
pendices. 



II. GENERALITIES AND FIRST PRINCIPLES 

Let the ideal (design) trajectory of a particle in an 
accelerator be a planar curve with curvature K. The 
Hamiltonian describing the motion of a particle in a nat- 
ural coordinate system attached to the orbit thus defined 
is is,]: 



(P, - qA^f - (P, - qA,f - q{l + Kx)As 



(1) 



where mpo is the rest mass of the particle. The guid- 
ing magnetic field can be represented as a gradient of a 
function 'il:(x,z;s) 

B = VV', (2) 

where the latter satisfies the Laplace equation 

V^t/- = 0. (3) 

Using the median symmetry of the machine, it is straight- 
forward to show that ip can be written in the form 



, , CL2X 

Tp ^ 1 ao + aix + + 



6o + bix + 
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^ + ico + ciX + ...)^ + . 



(4) 

Inserting the above expression into the Laplace equation 
^ , one readily finds relations between the coefficients bk 
and Cfc on one hand and on the other: 



^0 = Oq + Kai + 02, 



(5) 



bi = —2KaQ — K'gq + a" — K^ai + Ka2 + as, (6) 



+ 02 + 2K^ai - 2K'^a2 + Ka^ + 04, (7) 

co^b'^ + Kbi+b2. (8) 

Prime in the above expressions implies differentiation 
with respect to the longitudinal coordinate s. The co- 
efficients flfc have a very simple meaning: 



ao = {B, 



a2 



ai 



dBz 
dx 



dx"^ 



(9) 



In other words, this implies that, provided the vertical 
component Bz of the magnetic field and its derivatives 
with respect to the horizontal coordinate x are known in 
the median plane, one can in principle reconstruct the 
entire field chart. 

The vector potential A can be represented as 

Ar^ — —zF{x, z; s), Az = xF{x, z] s), A^ — G{x, z: s), 

(10) 

where the Poincare gauge condition 



xA^ -i- zAz = 0, 



(11) 



written in the natural coordinate system has been used. 
From Maxwell's equation 



&2 = &K^a'^ + GKK'a'Q - AKa'l ~ 2K'a[ 



B = V X A, 



(12) 
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we obtain 



2F + {xd^, + zd,)F = B„ 



(13) 



Kx — — 
j-^-j^G+{xd^ + zd,)G = zB,-xB,. (14) 

Applying Euler's theorem for homogeneous functions, we 
can write 



(15) 



[Co + Cix - 



4! 



(18) 



Inserting the above expansion into the Laplace equation 
for if, we obtain similar relations between Bk and Ck on 
one hand and Ak on the other, which are analogous to 
those relating b^, cj. and a^. 

We consider the canonical transformation, specified by 
the generating function 



S: 



z,T,P,,P„E;s] ^ xP., + zP, + TE 



1 Kx\ 

- H \B^^^ 

3 4 y " 



(16) 



where 



+ q j dTip{x,z,T; s). 



(19) 



(20) 



G 



zGx — xGz 
1 + Kx 



Here u — (a;, z) and si'^'' denotes homogeneous poly- 
nomials in X and z of order fc, representing the corre- 
sponding parts of the components of the magnetic field 
B = {Bx, Bz, Bs). Thus, having found the magnetic field 
represented by equation (l4|), it is straightforward to cal- 
culate the vector potential A. 

The accelerating field in AVF cyclotrons and FFAG 
machines can be represented by a scalar potential (the 
corresponding vector potential A = 0). Due to the me- 
dian symmetry, we have 



(17) 



is a canonical variable canonically conjugate to 7i. The 
relations between the new and the old variables are 

u^^^^^u, u^{x,z), r=^=r, (21) 
dPu 9E 



dS2 ^ 
ou 



J ATEu{ 



x,z,T; s) 



dip 



Pu-qEu{x,z,T-s), E^ = —^, (22) 



(/? = -I- Aix + 



- So + Six 
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'H = ^ = E + qip{x, z, T; s) = rripQ^c^ + qp{x, z, T; s). 

(23) 

The new Hamiltonian acquires now the form 



H^^{l + Kx)sl^~ ml^c^ - (Px ~ qEx - qA^)' - (p, - qE, - gA,) ' - q{l + Kx)(^As + E^j , (24) 



where 



/ 



Es = / dTE,{x,z,T-s) 



1 f ^^dp{x,z,T;s) 



1 + Kx 



ds 



We introduce the new scaled variables 



P, 



P P 



e = cT, 7 



E 



E 



Po "^poc Ep mpoc2 

(25) The new scaled Hamiltonian can be expressed as 



(26) 
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Pa 



— = -{l + Kx)Jj^-l~(^P,- qE, - gA,) - (P, - qE, - qA,^ - q{l + Kx) (a, + E^j , (27) 



r 



where 



Po 



(28) 



The quantities E^ and Ez can be neglected as compared 
to the components of the vector potential A, so that 



H = l3-i{\ + Kx) 



-^l-iP^-qA^ -iPz-qA,) -qA, 



where now 

- q q 



P PlPo ' 



P,, 



P. 



Pu 



P PlPo 



u = {x, z). 



(30) 

Since P„ and u are small deviations, we can expand the 
square root in power series in the canonical variables a;, 
Px and z, Pz- Tedious algebra yields 



H = H0+H1+H2 + H3 + H4 + (31) 



-q{l + Kx)E,, 



(29) 



Ho = -(3"f-q{l + Kx)E,, (32) 
Hi = p-fiqao - K)x, (33) 

H2 = ^ (pI + pI)+1 [{Kao + ai)x2 - a^z'] , (34) 



H3 = ^kx(pI + p: 



■{pI + pI) 



^-^{zP^-xPz) 



Kai 



Kai + a2 + y ]xz'^ 



(35) 



i?4 



p. 



— (Kao + 3a^){zP^~xPz) + — [x +z) 



( Ka2 03 , 4 



Ka2 



03 
3 



Kbn 



2 2 , 4 
- XZ +-Z 



(36) 



The Hamiltonian decomposition (|3ip represents the 
milestone of the synchrobetatron formalism. For in- 
stance, governs the longitudinal motion, Hi describes 
linear coupling between longitudinal and transverse de- 
grees of freedom and is the basic source of dispersion. 
The part H2 is responsible for linear betatron motion 
and chromaticity, while the remainder describes higher 
order contributions. 



III. THE SYNCHRO-BETATRON FORMALISM 
AND THE REFERENCE ORBIT 

In the present paper we consider a FFAG lattice with 
polygonal structure. To define and subsequently calcu- 
late the stationary reference orbit, it is convenient to use 
a global Cartesian coordinate system whose origin is lo- 
cated in the center of the polygon. To describe step by 
step the fraction of the reference orbit related to a partic- 
ular side of the polygon, we rotate each time the axes of 
the coordinate system by the polygon angle 6p = 2'k/Nl, 
where A''^ is the number of sides of the polygon. 

Let Xe and Pg denote the reference orbit and the ref- 
erence momentum, respectively. The vertical component 
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of the magnetic field in the median plane of a perfectly 
linear machine can be written as 



B,iXe;s)=aiis)[X,-X,-dis)l 



ao{Xe;s) = B^{Xe;s) 



(37) 



where s is the distance along the polygon side, and Xc is 
the distance of the side of the polygon from the center of 
the machine 



Xr = 



2tan(ep/2)' 



(38) 



Here Lp is the length of the polygon side which actually 
represents the periodicity parameter of the lattice. Usu- 
ally Xc is related to an arbitrary energy in the range from 
injection to extraction energy. In the case of EMMA it 
is related to the 15 MeV orbit. The quantity d{s) in 
equation ([57)1 is the relative offset of the magnetic cen- 
ter in the quadrupoles with respect to the corresponding 
side of the polygon. In what follows [see equations ((47|) 
and (|5D|) ] dp corresponds to the offset in the focusing 
quadrupoles and d^ corresponds to the one in the defo- 
cusing quadrupoles. Similarly, ap and Ud stand for the 
particular value of ai in the focusing and the defocusing 
quadrupoles, respectively. 

A design (reference) orbit corresponding to a local cur- 
vature K{Xf,; s) can be defined according to the relation 



K{X,;s) 



B,{X,-s) 



(39) 



where 7e is the energy of the reference particle. In terms 
of the reference orbit position Xc{s) the equation for the 
curvature can be written as 



(l + Xf)'/'s,(Xe;s), (40) 



where the prime implies differentiation with respect to s. 

To proceed further, we notice that equation (|40|) pa- 
rameterizing the local curvature can be derived from an 
equivalent Hamiltonian 

H,{X„P,-s) = -^^f^^ -q j AX,B,[X,-s) 



Taking into account Hamilton's equations of motion 
P. 



X'^ 



and using the relation 



Pr 



PeleX'c 



(43) 



we readily obtain equation (|40p . Note also that the 
Hamiltonian ([1T|) follows directly from the scaled Hamil- 
tonian ^ with a; = 0, P:r = -Pe, Pz = 0, = = 



and the accelerating cavities being switched off respec- 
tively. 

Hamilton's equations of motion (|42p can be linearized 
and subsequently solved approximately by assuming that 

Pe « I3ele. (44) 

Thus, assuming electrons (q = — e), we have 

eai(s) 



iX,^Xc-d{s)). 



(45) 

The three types of solutions to equations (|45p are as fol- 
lows: 

Drift Space 



Xc = Xq + (s - So), 

Pe7e 



Pe = Po, (46) 



where Xq and Pq are the initial position and reference 
momentum and s is the distance in longitudinal direction. 
Focusing Quadrupole 

Xe = Xc + dp + [Xq — Xc — dp) COSUJp{s ~ So) 



Pn 



sinujp{s — So), 



(47) 



Pe = -PeleUJpiXo - Xc - dp) sinwi^(s - So) 



where 



+ Po cosa;_F(s — sq), 



2 eap 



PoPele 

Defocusing Quadrupole 
Xc = Xc + do + {Xo - Xc - do) coshujois - so) 

D 

smhujp){s — So), 



(48) 
(49) 



(50) 



Pe = f3elei^D{Xa - Xc - do) sinhw£)(s - So) 



(41) 

P^ = gP.(Xe;s), (42) where 



Po cosh wd(s - So), 



eap 

PoPele 



(51) 



(52) 



In addition to the above, the coordinate transformation 
at the polygon bend when passing to the new rotated 
coordinate system needs to be specified. The latter can 
be written as 



Xp — Xc 



Xo — Xc 



cos Op - Po sin Qp/(3cie 
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Pe = Pele tail 



0„ + arctan ( ° 



(53) 



Once the reference trajectory has been found the cor- 
responding contributions to the total Hamiltonian (|3ip 
can be written as foUows 

iIo = -/37+^(^)/desin^(e), (54) 



(55) 



^2 = ^ (Pl + P,') +^ [(.g + MeK^)x' - g^] , (56) 



+ ^(2^3-3^?^), (57) 
6 



-"4 = :r-:—z , 



8r3 



24 



Pele 



(66) 



(67) 



The longitudinal part of the reference orbit can be iso- 
lated via a canonical transformation 

F2 (x,p^, z,p^, T, -q; = xp^ + zp^ + (r s) ^77 + , 

(68) 



(J = r -h s, 



(69) 



where a is the new longitudinal variable and 77 is the 
energy deviation with respect to the energy 7e of the 
reference particle. 



P2 + P 



-z . 



SP^-f^ 24 
Here, we have introduced the following notation 
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qai 
Po ' 



(58) 



(59) 



Moreover, Z is the charge state of the accelerated parti- 
cle, A is the mass ratio with respect to the proton mass 
in the case of ions, and ^(O) is the phase of the RF. 
For a lepton accelerator hke EMMA, A = Z = 1. In 
addition, (dAE/ds) is the energy gain per unit longitu- 
dinal distance s, which in thin lens approximation scales 
as AE/As, where As is the length of the cavity. It is 
convenient to pass to new scaled variables as follows 



Pu 



Pu 



h = 



7 



(60) 



IV. DISPERSION AND BETATRON MOTION 

The (linear and higher order) dispersion can be intro- 
duced via a canonical transformation aimed at canceling 
the first order Hamiltonian Hi in all orders of r]. The 
explicit form of the generating function is 

G2{x,px,z,pz,a, fj; s) = arj + 1% + xp^ 



Y,ri''[^Ms)-P.ms)+Sk{s)], (70) 

k=l 



X ^x + ^ff'Vk, 
fe=i 



= P. + J2^'''^k, (71) 



fe=i 



r^p.e. r..|2_.^ftW-^^ (01) 

Thus, expressions (f54|) - ((58|) become 

Z /'dAE\ f 
i/o^-re + ^(^jydrsin0(.), (62) 

Hi = -(Fe - l)Kx, (63) 

H2 = ^^{pi+Pl) + 1 ^^'^ + - '^^^ ' 

Hs-^{pi+pi) + ^ {2x' - ix^) , (65) 



CSO 

<^ = ^ + Yl kt'^'PkPx - XkX) 
fe=l 



- J2 ^t'^ [Sk+XkY. V"'Vra ■ (72) 
k=l \ m=l J 

Equating terms of the form xrf^ and Px'^T' in the new 
transformed Hamiltonian, we determine order by order 
the conventional (first order) and higher order disper- 
sions. The first order in 77 (terms proportional to xrj and 
PxTf) yields the well-known result 

V[ = Xi, X[ -f (ge + K'')Vi = K. (73) 

Since in the case of vanishing betatron motion 
(x = 0, Px = 0) the new longitudinal coordinate t? 



7 



should not depend on the new longitudinal canonical con- 
jugate variable 77, the second sum in equation ([7^ must 
be identically zero. We readily obtain iSi = 0, and 



S2 



XiVi 



In second order we have 



K 



and in addition the function 53(3) is expressed as 

S3 = -liXiV2 + 2X2ri). 



(74) 

(75) 
(76) 

(77) 



Close inspection of equations ([75]) , ((75l) and ([75)1 shows 
that Vi is the well-known linear dispersion function, V2 
stands for a second order dispersion and so on. Up to 
third order in rj the new Hamiltonian describing the longi- 
tudinal motion and the linear transverse motion acquires 
the form 

^ iC^fj' , iC2v' , Z fdAE\ 

(78) 
(79) 



H2 = ^{pI+pI) 



where 



KVi 



XI 



For the sake of generality, let us consider a Hamiltonian 
of the type 



u—{x,z) 



(81) 



A generic Hamiltonian of the type ([5T|) can be trans- 
formed to the normal form 



u—{x^z) 



(82) 



by means of a canonical transformation specified by the 
generating function 

.f2(x,P„z,P.;,s)= E f^-^y (83) 



u—{x^z') 



Here the prime implies differentiation with respect to the 
longitudinal variable s. The old and the new canonical 
variables are related through the expressions 



1 



Pu 



(Pu-auU). (84) 



14 16 

Energy [MeV] 



FIG. 1: Horizontal betatron tune for the EMMA ring as a 
function of energy. 



The phase advance Xu{s) and the generalized Twiss pa- 
rameters au(s), Puis) and ju{s) are defined as 



da,. 



ds 



Pu 



as 



ds 



(85) 



(86) 



(87) 



The third Twiss parameter 7u(s) is introduced via the 
well-known expression 



Pulu - Q!„ = 1. 



The corresponding betatron tunes are determined ac- 
cording to the expression 



2^ J i3u{e) 



(89) 



Typical dependence of the horizontal and vertical be- 
tatron tunes on energy in the EMMA non scaling FFAG 
is shown in Figures [T] and [2] 



V. ACCELERATION IN A NON SCALING 
FFAG ACCELERATOR 

The process of acceleration in a non scaling FFAG ac- 
celerator can be studied by solving Hamilton's equations 
of motion for the longitudinal degree of freedom. The lat- 
ter are obtained from the Hamiltonian (|4ip supplemented 
by an additional term [similar to that in equation (|54p ]. 
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14 16 
Energy [MeV] 



FIG. 2; Vertical betatron tune for the EMMA ring as a func- 
tion of energy. 




10 12 14 16 18 20 

Energy [MeV] 

FIG. 3: Time of flight as a function of energy for a single 
0.394481 meter EMMA ceU. 



which takes into account the electric field of the RF cav- 
ities. They read as 



de 

d^ 



7 



^^2^2 _ p2 ' 



d7 
dl 



ZelJc \rA 



Sfc)sin 



fc=i 



(90) 



(91) 



Here Uc is the cavity voltage, ujc is the RF frequency, 
Nc is the number of cavities and ipk is the corresponding 
cavity phase. 

One could use the results obtained in the previous sec- 
tion with the additional requirement that the phase slip 
coefficient /Ci averaged over one period vanishes. Instead, 
we shall use an equivalent but more illustrative approach. 
The path length in a FFAG arc and therefore the time of 
flight O is often well approximated as a quadratic func- 
tion of energy. The acceleration process is then described 
by a longitudinal Hamiltonian, which contains terms pro- 
portional to the zero-order (conventional phase slip) fac- 
tor and first-order phase slip factor. It usually suffices 
to take into account only terms to second order in the 
energy deviation 



(92) 



as suggested by Figure [31 

Here 7™ corresponds to the reference energy with a 
minimum time of flight. Provided the time of flight Qi at 
injection energy 7^ and the time of flight Q„i at reference 
energy 7^ are known, the constants entering equation 
([^ can be expressed as 



A = 



e,. 



(7m - l^f 



60 = e„, - ^7,2 



(93) 



Next, we pass to a new variable 

7 = 7- 7m, 6 = 0r; 



(94) 



similar to the variable rj introduced in the previous sec- 
tion. Then, Hamilton's equation of motion ((90)) can be 
rewritten in an equivalent form 



de 



^7' 



ds 



(95) 



In what follows, it is convenient to introduce a new 
phase ip and the azimuthal angle 9 along the machine 
circumference as an independent variable according to 
the relations 



ds = Rde, 



NlLp 
2ti 



(96) 



It is straightforward to verify (see the averaging proce- 
dure below) that the necessary condition to have accel- 
eration is 



UJcNl\Q^\ 

2ttc 



(97) 



where h is an integer (a harmonic number). Averaging 
Hamilton's equations of motion 



dip 



haj , 



d7 
d^ 



ZeUc 



A 



9k) sin {ip- ipk), 



fe=i 



we rewrite them in a simpler form as 



dp 
'dB 



harj 



d7 
d^ 



A sin (/?, 



(98) 



(99) 



(100) 



9 



(108) 



In the above expressions 2Fi(c(, denotes the 

Gauss hypergeometric function of the argument x. This 
case is illustrated in Figure [H 

In the general case where Hq ^ 0, we have 

by/Qic \6 2 2 6 ai c J b 

where 



FIG. 4: An example of the so-called serpentine acceleration 
for the EMMA ring for the central trajectory, where the lon- 
gitudinal Ho = 0. The harmonic number is assumed to be 11, 
with the RF wavelength 0.405m. The parameter a from Eq. 
is taken to be 2.686310"^ 



Ho 
X ' 



■.Fi 



11 1. 4 J3 

3' 2' 2' 3' 



Ho 
A ' 



-ll. 
c 



(110) 



(111) 



where 



(fi = -(f - he + tpo, 



A = ^ , 101 



Here now, F\(a\ /3, 7; 5; x, y) denotes the Appell hyperge- 
ometric function of the arguments x and y. The phase 
portrait corresponding to the general case for a variety of 
values of the longitudinal Hamiltonian is illustrated 
in Figure [5l 



%po = arctan ( ^ ) , (102) 



VI. CONCLUDING REMARKS 



Ac^^ cos {hOk + (^fe), As = ^ sin {hOu + (/?&). 



fe=i 



(103) 

The effective longitudinal Hamiltonian, which governs 
the equations of motion pOO|l can be written as 



Hq = —Y + Acosv?. 



(104) 



Since the Hamiltonian (|104p is a constant of motion, the 
second Hamilton equation (|100p can be written as 



d-y 
d^ 



±AWl 



A2 



Ho 



-T 



(105) 



Let us first consider the case of the central trajectory, 
where Ho — 0. It is of utmost importance for the so 
called gutter acceleration. Equation pOSp can be solved 
in a straightforward manner to give 



-J pf^ 1.7. C 



where 



(106) 



(107) 



Based on the Hamiltonian formalism, the synchro- 
betatron approach for the description of the dynamics of 
particles in non scaling FFAG machines has been devel- 
oped. Its starting point is the specification of the static 
reference (closed) orbit for a fixed energy as a solution of 
the equations of motion in the machine reference frame. 
The problem of dynamical stability and acceleration is 
sequentially studied in the natural coordinate system as- 
sociated with the reference orbit thus determined. 

It has been further shown that the dependence of the 
path length on the energy deviation can be described in 
terms of higher order (nonlinear) dispersion functions. 
The method provides a systematic tool to determine the 
dispersion functions to every desired order, and repre- 
sents a natural definition through constitutive equations 
for the resulting Twiss parameters. 

The formulation thus developed has been applied to 
the electron FFAG machine EMMA. The transverse and 
longitudinal dynamics are explored and an initial at- 
tempt is made at understanding the limits of longitudinal 
stability of such a machine. 

Unlike the conventional synchronous acceleration, the 
acceleration process in FFAG accelerators is an asyn- 
chronous one in which the reference particle performs 
nonlinear oscillations around the crest of the RF wave- 
form. To the best of our knowledge, it is the first time 
that such a fully analytic theory describing the accelera- 
tion in non scaling FFAGs has been developed. 
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FIG. 5: Examples of serpentine acceleration for the EMMA 
ring, with varying value of the longitudinal Hamiltonian. The 
limits of stability are given at values of the longitudinal Hamil- 
tonian of ±0.31272, corresponding to either a phase at 
lOMeV, or a TT phase at 20MeV. 



APPENDIX A: CALCULATION OF THE 
REFERENCE ORBIT 

The explicit solutions of the linearized equations of mo- 
tion ((45|) can be used to calculate approximately the ref- 
erence orbit. To do so, we introduce a state vector 




(Al) 



The effect of each lattice element can be represented in 
a simple form as 



(A2) 



Here Zi„ is the initial value of the state vector, while 
Zo«t is its final value at the exit of the corresponding 
element. The transfer matrix M.ei and the shift vector 
Agj for various lattice elements are given as follows: 
1. Polygon Bend. 

Within the approximation (|44p considered here we can 
linearize the second of equations (|53p and write 



l/cos8p — Xctan6p/(/3e7eCos9p) 
l/cos^Op 

Xc{l - 1/ cos Op) 
/3e7e tan 9p 



(A3) 



2. Drift Space. 



Mo = 



1 Lo/f3ele 







Ao = 0, (A4) 



where Lq is the length of the drift. Every cell of the 
EMMA lattice includes a short drift of length Lq and a 
long one of length Li. 
3. Focusing Quadrupole. 

The transfer matrix can be written in a straightforward 
manner as 

Mf=\ I , (A5) 

-/3e7e(jji? sin (wi^ii?) cos{ujfLf) 



[Xc + (i_F)[l — COS {uofLf)] 
PeleUJF{Xc + dF) sin (lofLf) 



(A6) 



where Lf is the length of the focusing quadrupole. 
4. Defocusing Quadrupole. 

The transfer matrix in this case can be written in anal- 
ogy to the above one as 



Md = 



An = 



C0sh(w£,LD) Smh{LUDLD)/iPelet^D) 
PeJei^DSmh{ujDLD) COsh(wDi_D) 

{Xc + dD)[l- cosh (ujdLd)] 
-Pele^oiXc + do) sinh {ujdLd) 



(A7) 
(A8) 



where Ld is the length of the defocusing quadrupole. 

Since the reference orbit must be a periodic function 
of s with period Lp, it clearly satisfies the condition 

Zout = Zi„ = Ze. (A9) 

Thus, the equation for determining the reference orbit 
becomes 

Z, ^ MZc + A, or Ze = (l - M) ^ A. (AlO) 



Here M. and A are the transfer matrix and the shift 
vector for one period, respectively. The inverse of the 
matrix \ — M. can be expressed as 



l-M 



cos'^ 9, 



1 + 1 - SpTW cos3 Op 



I-M22 M12 
M21 l-Mu 



(All) 



For the EMMA lattice in particular, the components 
of the one period transfer matrix and shift vector can be 
written explicitly as 
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CpCD + [ — - LqLiOJfUJD ] SFSD + {Lo + LijOJDCFSD — LiLOFSFCD 
U)F 



(A12) 



Ml2 



1 r/Lo + ^i 



V y. \ p 



Xctpj Cfcd + 



LqLiUJfOJD —\Xctp — — - 

LUf J UJDCp 



SfSd 



UJDCp 



— (Lo + Li)uj£iXctp 



, / 1 , , ^ , LoLiujf\ 

CfSd + \- LiLOF^ctp SfCd 

\UJFCp Cp J 



(A13) 



M21 = —{ljfsfcd + LqljfujdSfsd - ojdCfsd), 



(A14) 



M22 = 



cfCd 



U)f \ ( io 1 

LoUJFUJDXctp SFSD + UJF Xctp SfCd — UlDXctpCFSD 

UJDCp J \ C, 



(A15) 



Ai ^ Xc + dp + [d-D — dF){cF — LiUfSf) + [ — - + dD 

Cn 



ujdsfsd 

LilofSfCd - cfCd - [Lo + Li)ujdCfSd 1- LqLiluflodSfSd 

UlF 



cfsd , SfCd Liujfsfsd j- ^ 

(Lq + LijCFCD ^ ! LqLiUJfSfCd 

ujd ujf ujd 



(A16) 



A2 = -PeleUJFidD - d_F)sF + Pe^e ( — + ) [uJfSfCD + UJfUJdLqSfSd - UJdCfSd) 



, « W UJFSFSD J- \ 

+ Peletp CfCd UJfLqSfCD ■ 

\ UJD J 



(A17) 



For the sake of brevity, the foUowing notations 

Cp = cos6p, CF — cos{ujfLf), cd = cosh (ujd Ld), 

(A18) 

<:p = tan0p, sf = sin (wf-^f), sd = sinh (wdL^)), 

(A19) 

have been introduced in the final expressions for the com- 
ponents of the one period transfer matrix and shift vec- 
tor. 



APPENDIX B: PHASE STABILITY IN FFAGS 

To study the stabihty of the serpentine acceleration in 
FFAG accelerators, we write the longitudinal Hamilto- 



nian (|104p in an equivalent form 

Ho = + cosip). (Bl) 
Hamilton's equations of motion can be written as 

Let >fai&) and Ja(0) be the exact solution of equations 
(|B2[) described already in Section V. Let us further de- 
note by ipi and Ji a small deviation about this solution 
such that ip — + and J = Ja + Ji- Then, the 
linearized equations of motion governing the evolution of 
if I and J I are 

^^=66JaJl, ^ = 5<piC0S(^a- (B3) 
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FIG. 6: Phase stability of the standard EMMA ring, for the 
central trajectory at Ho ~ 0. The errors are given as O.lMeV 
in energy and 1.3° in phase. 



The latter should be solved provided the constraint 

3Jf Ji-(^isin(p, = 0, (B4) 



following from the Hamiltonian jBTj) holds. Differentiat- 
ing the second of equations (|B3p with respect to 9 and 
eliminating tpi, we obtain 



dVi Gb^Ho 



JaJi + Ibb^J^Ji - 



de^ 



A 



(B5) 



Next, we examine the case of separatrix acceleration with 
Hq = 0. In Section V we showed that to a good ac- 
curacy the energy gain [Ja{S) = bd + Ji] is linear in the 
azimuthal variable 9. Therefore, equation (|B5|) can be 
written as 



The latter possesses a simple solution of the form 



(B6) 



Ji 



\Ja\ 




C2yi 




(B7) 



where Ja{z) and 3^q(z) stand for the Bessel functions of 
the first and second kind, respectively. In addition the 
constants C\ and C2 should be determined taking into 
account the initial conditions 



dJi(J.) 

dJa 



= (pi(Ji) COS(/?i 



J\{Ji) — Jli 



(B8) 
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